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ABSTRACT 



We carried out global three-dimensional magneto-hydrodynamic simulations 
of dynamo activities in galactic gaseous disks without assuming equatorial sym- 
metry. Numerical results indicate the growth of azimuthal magnetic fields non- 
symmetric to the equatorial plane. As magneto-rotational instability (MRI) 
grows, the mean strength of magnetic fields is amplified until the magnetic pres- 
sure becomes as large as 10% of the gas pressure. When the local plasma f3 
(— Pgas/Pmag) becomes less than 5 near the disk surface, magnetic flux escapes 
from the disk by Parker instability within one rotation period of the disk. The 
buoyant escape of coherent magnetic fields drives dynamo activities by gener- 
ating disk magnetic fields with opposite polarity to satisfy the magnetic flux 
conservation. The flotation of the azimuthal magnetic flux from the disk and 
the subsequent amplification of disk magnetic field by MRI drive quasi-periodic 
reversal of azimuthal magnetic fields in timescale of 10 rotation period. Since 
the rotation speed decreases with radius, the interval between the reversal of 
azimuthal magnetic fields increases with radius. The rotation measure computed 
from the numerical results shows symmetry corresponding to a dipole field. 

Subject headings: Galaxy: disk - Magnetohydrodynamics(MHD) - Dynamo 
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Introduction 



The origin and nature of magnetic fields in spiral galaxies is an open question. 
Magnetic fields in spiral galaxies have been int ensively studied through observations 



of linear po 



Beck et al. 



arization of synchrotron radiation ( ISofue. Fujimoto. fc Wielebinski 



19861 : 



19961 ) . The typical magnetic field strength is a few //G in the Milky Way near 



the Sun, and in spiral galaxies. T 



to unity ( iRand fc Kulkarni 



r e mean ratio of gas p ressure to magnetic pressure is close 



1989 



Crocker et al. 



201(1 ). suggesting that magnetic fields 



in spiral galaxies are strong enough to affect the motion of the interstellar matter in the 
gaseous disk. 

Since the galactic gas disk rotates differentially, the magnetic field is amplified faster 
near the galactic center than in the outer disk, and the magnetic field strength can exceed 
O.lmG in the galactic center. An observational evidence of such strong magnetic fields has 
been obtained in the radio arc near the galactic center, where the field strength of about 
lmG is estimated from synchrotron radiation. Ma ny threads perpendicular to the galacti c 



plane are obs erved indicating strong vertical field ( lYusef-Zadah. Morris. &: Chance 



1984 



Morris 



19901 ) . On the other hand, observations of near infrared p olarization su 



that horizontal fields of 0.1 — lmG are dominant inside the disk (INovak et al. 



ggested 



2000; 



Nishiyama et al. 



20091). 



Taylor. Stil. fc Sunstrum 



( 120091 ) presented results of all sky distribution of rotation 
measure (RM), which indicate that the magnetic fields along line-of-sights are roughly 
consistent with the dipole-like models in which the azimuthal magnetic fields are anti- 
symmetric with respect to the equatorial plane, rather than quadrupole-like models in 
which the azimuthal magnetic fields show symmetry with respect to the equatorial plane. 
They also showed that the RM distribution has s mall scale variation s, which would be 



linked to turbulent structure of the gaseous disk. 



Han et al. 



( 120021 ) showed a view of 
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RM distributions from the North Galactic pole using the results of pulsar polarization 
measu rements. They indicated that the mag netic field direction changes between the spiral 



arms (IHan et al. 



2002 



Han fc Zhang 



20071 ) 



The magnetic fields have to be maintained for 10 £ 



10 years, be cause most of 



spiral galaxies have magnetic fields with similar strength ( jBeck 



20121 ). 



Parker 



(1971) 



considered the self-sustaining mechanism of magnetic fields called the dynamo activity. In 
the framework o f the kinematic dynamo, where the back-reaction by enh anced magnetic 
fields is ignored, 



Fujimoto fc Sawa 



(119871 ) and ISawa fc Fuiimoto I (|1987l) obtained the 



Balbus fc Hawley 



(1991) 



global structure of galactic magnetic fields. On the other hand, 
pointed out importance of the magneto-rotational instability (MRI) in differentially rotating 
system, and showed that the back-reaction of the enhanced magnetic fields cannot be 



ignored. Since ga" 



actic gaseous disk is a d ifferentially rotating system, magneto-rotational 



instability grows (IBalbus fc Hawley 



19911 ) 



In this paper, we present the results of global three-dimensional magneto-hydrodynamic 
simulations of galactic gas disks assuming an axisymmetric gravitational potential and 
study the dynamo activities of the galactic gas disk. 

In order to take into account the back reaction of magnetic fields on the dynamics of 
gas disk, we have to solve the magneto-hydrodynamic equations in which the induction 
equation is coupled with the equation of motion. This approach is called the dynamical 
dynamo. Nishikori, Machida & Matsumoto (2006) presented results of three dimensional 
magneto-hydrodynamic (MHD) simulations of galactic gaseous disks. They showed that, 
even when the initial magnetic field is weak, the dynamical dynamo amplifies the magnetic 
fields and can attain the current field strength observed in galaxies. They also showed that 
the direction of mean azimuthal magnetic fields reverses quasi-periodically, and pointed 
out that reversal of magnetic field is driven by the coupling of MRI and Parker instability. 
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When the magnetic field is amplified by MRI, the time scale of buoyant escape of magnetic 
flux due to Parker instability becomes comparable to the growth time of MRI, and limits 
the strength of magnetic fields inside the disk. 



Since 



Nishikori et al. 



( 120061 ) assumed symmetric boundary condition at the equatorial 



plane of the disk, they could not simulate the growth of anti-s ymmetric mode of Pa rker 



1988|). The 



instability, in which magnetic fields cross the equatorial plane (IHoriuchi et al. 
present paper is aimed at obtaining a realistic, quantitative model of the global magnetic 
field of the Milky Way Galaxy based on numerical 3D MHD simulations of magnetized 
galactic gaseous disk. In section 2, we show the initial setting of the simulation. We present 
the results of numerical simulations in section 3. The result will be used to compute the 
distribution of RM on the whole sky in section 4. 



2. Numerical Models 



2.1. Basic equations 



We solved the resistive MHD equations by using modified Lax-Wendroff scheme with 
artificial viscosity (Rubin & Burstein 1967, Richtmyer & Morton 1967). 



dp 
dt 



+ V ■ pv = 



dv 11 

— + (v ■ V)u = — VP - V0 + (V x B)xB 

at p 4iip 



d 



£? 2 



di\2 pV + 8^ + 7 ^l 



-pv 2 + — — v + —E x B 
2 7- 1/ 47r 



dB 

~dt 



V x (v x B - r]V x B) 



(1) 
(2) 

-pvV(p (3) 
(4) 



where p, P, v, B, and 7 are the density, gas pressure, velocity, magnetic field, the specific 
heat ratio, respectively, and is gravitational potential proposed by Miyamoto & Nagai 
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(1975). The electric field E is related to magnetic field B by Ohm's law. We assume the 
anomolous resistivity 77 

{r]o(v d /v c -l) 2 , v d >v c , 
(5 
0, v d < v c , 

which sets in when the electron-ion drift speed v d = j/p, where j is the current density, 

exceeds the critical velocity v c . 



2.2. Simulation Setup 



We assumed an equilibrium to rus threaded by weak toroidal magnetic fields 



(jOkada. Fukue fc Matsumoto 



1989). The torus is embedded in a hot, non- rotating 
spherical corona. We assume that the torus has constant specific angular momentum L 
at w = cub and assume polytropic equation of state P = Kp 1 , where K is co nstant. The 



density distribution is same as that for equation (8) in 



Nishikori et al. 



(j2006h 



We consider a gas disk composed of one-component interstellar gas. Radiative cooling 
and self-gravity of gas are ignored in this paper. 

We adopt a cylindrical coordinate system (w,(p,z). The units of length and velocity 



in this paper are zuq = 1 kpc and Vq = (GM/wq 



vl/2 



207km/s, respectively, where 



M = 10 Mq. The unit time is to = ^o/ v o- Other numerical units are listed in Table 1. 

The numbers of grid points we used are (iV w , N v , N z ) = (250, 128, 640). The grid size 
is Avu/wq = 0.05, Az/mo = 0.01 for < vj/vjq < 6.0, \z\/wq < 2.0, and otherwise they 
increase with w and z, respectively. The grid size in azimuthal direction is Aip = 27r/128. 
The outer boundaries at vd — 56g7o and \z\ = 10wo are free boundaries where waves can 
be transmitted. We applied periodic boundary conditions in the azimuthal direction. An 
inner absorbing boundary condition is imposed at r = a/gj 2 + z 2 = r in = 0.8cj7 , since the 
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gas accreted to the central region of the galaxy will be converted to stars or swallowed 
by the black hole. To illuminate the effect of boundary condition at the equatorial 
plane, we calculate two models: one is model PSYM which imposes sy mmetric boundar y 
condit ion at the equatorial plane. This model is identical with that in 



Nishikori et al. 



( 120061 ). The other is model ASYM, in which we include the equatorial plane inside the 



simulation region, and imposed no boundary condition at the equatorial plane. In this 
paper, we adopted model parameters = IOzuq, 0h = 100 at (va,z) = (zu\>,0), 7 = 5/3, 
v c = lOOfo, and rjo = O.Ivqiuo. The thermal energy of the torus is parameterized by 
Eth = c lh/(^ v o) = Kpl~ 1 / v o = 0-^5 where c s b and pb are the sound speed and density 
at (vj,z) = (t*7b,0), respectively. The initial temperature of the torus is T ~ 2 x 10 5 K. 
This mildly hot plasma mimics the mixture of hot (T ~ 10 6 K) and warm (T ~ 10 4 K) 
components of the interstellar matter, which occupy large volume of the galactic disk. 



3. Numerical results 

3.1. Time evolution of the galactic gaseous disk 

Figured^, shows the time evolution of the mass accretion rate measured at w/wq = 2.5. 
After t/to ~ 200, the gas starts to inflow and the mass accretion rate increases linearly. 
Subsequently (t/t > 600) the accretion rate saturates and becomes roughly constant. The 
evolution of the mass accretion rate in model ASYM (black) is similar to that in model 
PSYM (gray), where the mass accretion rate for model PSYM is doubled since the model 
PSYM includes only the region above the equatorial plane. 

Figures [lb, [lb and [Tbl show the time evolution of the plasma (3 and magnetic energy of 
each component averaged in the region where 2 < w/wq < 5, \z\/vjq < 1, and < (p < 2tt, 
respectively. It is clear that plasma (3 decreases in the linear stage as the magnetic energy 
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increases. In both models, the averaged magnetic energy first increases exponentially. After 
that, it saturates and becomes roughly constant. Time evolution of magnetic fields is similar 
between model ASYM and model PSYM. Time evolutions of the growth and saturation in 
the mass accretion rate are quite similar to those of the averaged magnetic energy. It means 
that the magnetic turbulence driven by MRI is the main cause of the angular momentum 
transport which drives mass accretion. In the remaining part of this section, we discuss the 
results of model ASYM. 

In order to check the magnetic fields structure, we analyze the time evolution of 



mean azimutha 



Nishikori et al. 



magn etic fields. The mean fields are computed by the same method as 



(120061 ) . Figure shows the time evolution of mean azimuthal magnetic 
fields. Black curve shows B v averaged in the region where 5 < cu/t^o < 6, < y? < 27r, 
and < z/vjq < 1, and gray curve shows B v averaged in the region where 5 < -to/wo < 6, 
< if < 27r, and 1 < z/wq < 3. Black and gray curves correspond to the disk region 
and halo region, respectively. The azimuthal magnetic fields reverse their direction with 
timescale t ~ 300to ~ 1.5Gyr. This timescale is comparable to that of the buoyant rise of 
azimuthal magnetic fields. The halo fields also reverse their directions with the same time 
scale as the disk component. 

The radial distribution of the mean azimuthal magnetic fields at t = lOOOto are shown 
in figure [2b. Black and gray curves show the averages in the disk (0 < ip < 2ir and 
< zjvjQ < 1), and in the halo (0 < (p < 2ir and 1 < z/w < 3), respectively. Since the 
magnetic fields in the disk becomes turbulent, the direction of azimuthal magnetic fields 
changes frequently near the equatorial plane. On the other hand, long wavelength magnetic 
loops are formed in the halo region. 

Figure [3^ shows the radial distribution of the specific angular momentum at equatorial 
plane. Although the rotation is assumed only inside initial torus (gray curve in figured), 
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torus deforms its shape from a torus to a disk by redistributing angular momentum (black 
curve in figure Eh)- Since the outflow created by the MHD driven dynamo transports 
angular momentum toward the halo, the halo begins to rotate differentially, as can be seen 
in iso-contours of the azimuthally averaged rotation speed (figure [3b). Figure [3^ and |3b 
show that the azimuthal velocity v v is almost constant and close to vq in the equatorial 
plane. The rotation period at zu/w = a is P rot ~ 27ia(zu /v ) = 2nat . 



3.2. Evolution of the magnetic field structure 

We show the 3D structure of mean magnetic fields in figure HI The colored surfaces 
represent iso-density surfaces, and the curves show magnetic field lines. Light brown curves 
show field lines passing through the plane at z/w$ = 5 and wjw^ < 1. The vertical 
magnetic fields around the rotation axis are produced by the magnetic pressure driven 
outflow near the galactic center. Color on the curves depicts the direction of azimuthal 
magnetic fields. Red curves show positive direction of azimuthal magnetic fields (counter 
clockwise direction) and blue shows negative direction (clockwise). The figure indicates 
that the azimuthal field reverses the direction with the radius. 



Figure H] indicates that long wavelength magnetic loops are formed around disk-halo 
interface. Therefore, we discuss the formation mechanism of the magnetic loops. Horizontal 
magnetic fields embedded in gravitationally stratified atmosphere become unstable against 
long wavelength undular perturbations. This undular mode of the magnetic buoyancy 



insta bility, called Parker instability, creates buoyantly rising magnetic loops ((Parker. 



19661 ) . Nonlinear growth of Parker instabil ity in gravitationa 



by magneto-hydrodynamic simulations by 



Matsumoto et al- 



l y stra tified disks was studied 



(11988h . They showed that 



magnetic loops continue to rise when /3 < 5. In we akly magnetized region w here /3 > 5 



Parker instability only drives nonlinear oscillations (IMatsumoto et al. 



1990). 
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In differentially rotating disks, MRI couples with Parker instability. The growth 
time scale of the non-axisymmetric MRI is ^mri ~ l/(0.1fl) ~ 10H/c s where c s is 
the sound speed. On the other hand, the growth time scale of Parker instability is 
tpi ~ 5H/va ~ 5^//3 , y/2(H/c s ) where va and 7 are the Alfven speed and the specific 
heat ratio, respectively. This time scale becomes comparable to the growth time of MRI 
when /3 ~ 5. Since magnetic turbulence produced by MRI limits the horizontal length 
of coherent magnetic fields, the growth of Parker instability can be suppressed in weakly 
magnetized region inside the disk. On the other hand, in regions where (3 < 5, nonlinear 
growth of Parker instability creates magnetic loops buoyantly rising to the disk corona 



(IMachida. Hayashi. fc Matsumoto 



2000 



Machida. et al. 



20091 ). 



Figure [5^ shows an example of a magnetic loo p at t = lOOOtn- Magne tic loops are 



identified by the same algorithm as that reported in 



Machida. et al. 



( 120091 ). Figure Eb, Eb, 



and [5H show the distribution of the vertical velocity, plasma (3, and density, respectively 
along the magnetic field line depicted in figure EK- A black circle shows the starting point 
of the integration of a magnetic field line. The positive vertical velocity indicates that the 
magnetic loop is rising. The value of plasma (3 decreases from the foot-point of the loop 
where (3 ~ 100 to the loop top where (3 < 5 (figure |5H) . The density also decreases from the 
foot-points to the loop top. These are typical structure of buoyantly rising magnetic loops 
formed by Parker instability. 

From the numerical results presented so far, the initial magnetic field was found to 
be amplified and transported to the halo region after several rotation periods. In order 
to investigate the time evolution of the azimuthal magnetic field, we show a butterfly 
diagram of the azimuthal magnetic fields at w/wq = 2 (figure [6^) • The horizontal axis 
shows time t/to, and the vertical axis shows the height z/wq. Color denotes the direction 
of the azimuthal magnetic fields. White curve shows the iso-contour where (3 = 5 above the 
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equatorial plane. We should note that there is no magnetic flux at w/vjq = 2 initially. As 
the time elapses, the accretion of the interstellar gas driven by MRI transports magnetic 
fields from vj/vjq = 10 to w/zo = 2. The amplified azimuthal magnetic fields change their 
directions quasi-periodically. Turbulent magnetic fields are dominant around the equatorial 
plane. 

When the magnetic tension becomes comparable to that of turbulent motion in the 
disk surface and coherent length of the magnetic field line increases, the magnetic flux 
buoyantly escapes from the disk due to Parker instability. We find that the instability 
significantly grows, when the plasma /3 decreases to around 5. The averaged escape velocity 
of the flux is about 25 km/s, approximately equal to the Alfven velocity. 

The periodic reversals of azimuthal magnetic fields can also be seen in correlation 
between the fields above and below the equatorial plane (figure EJ^b)). White color shows 
the region where B v (z) ■ B^—z) > 0, and black shows the region where B v (z) ■ B^{—z) < 0. 
We can see fine structures of turbulence near the equatorial plane. In the halo region where 
z/wq > 1, the color of outgoing flux changes quasi-periodically. It means that the disk 
magnetic topology changes between symmetric state and anti-symmetric state. 

Finally, we check the spatial distribution of the azimuthal magnetic fields. The 
azimuthal magnetic fields in w — z plane averaged in the region < (p < 2tt at t/t = lOOOis 
shown in figure [7k.. We can see that a few /zG fields are distributed as high as at z/w = 5. 
Figure [7b, and [7b show the distribution of the azimuthal magnetic fields in w — <p plane, 
where the azimuthal magnetic fields are averaged in the halo (1 < z/vjq < 1.5) and in the 
disk (0 < z/vjq < 0.5), respectively. Since turbulent magnetic fields are dominant inside the 
disk, short magnetic filaments are formed (figure [7b ). Long and strong azimuthal magnetic 
sectors are formed in the halo region, because the plasma in the halo region becomes 
lower than that in the disk due to the density decrease (figure [7b). Therefore, magnetic 
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tension of buoyantly rising magnetic loop exceeds the ram pressure by turbulent motion 
around the disk surface and in the halo. These magnetic loops buoyantly escape from the 
disk by Parker instability and create the rising magnetic fluxes in the butterfly diagram 
shown in figure [6k. 



3.3. Dependence on the azimuthal resolution 

In order to study the dependence of numerical results on azimuthal resolution, we 
carried out simulations in which we used 64 (model ASYM64) and 256 (model ASYM256) 
grid points in < ip < 2tt. Other parameters are the same as those in model ASYM. 

Figure [Hk shows the time evolution of the magnetic energy averaged in the same region 
as figure [lb. Black, dashed and gray curves show model ASYM, ASYM256, and ASYM64, 
respectively. This panel indicates that the saturation level of the magnetic energy near the 
equatorial plane slightly decreases as the resolution increases because shorter wavelength 
turbulence inside the disk can be resolved better in the model with higher azimuthal 
resolution, so that more magnetic energy dissipates in the nonlinear stage. The time 
evolution of the azimuthal magnetic fields averaged in the region where 5 < w/vj$ < 6, 
< z/wq < 1, and < (p < 2ir is shown in figure [Hb. Colors are same as figure [Hh. 
The amplitude of the azimuthal magnetic fields and the time scale of the field reversal 
are almost similar in model ASYM and ASYM256. In order to resolve the most unstable 
wavelengt h of MRI, we need more than 20 grid points per disk thickness for azimuthal 



direction (lHawley. Guan. fc Krolik 



20111 ) . Although ASYM does not satisfy this condition, 
the qualitative behavior of the results for ASYM is identical to that of ASYM256, which 
satisfies the condition (see figure Eb)- For example, the ratio of the radial magnetic energy 
to the azimuthal magnetic energy (B^/B^) is 0.12 in model ASYM and 0.13 in model 
ASYM256. Hence, we conclude that numerical results of model PSYM and ASYM in 



13 



this paper obtained by simulations with 128 grid points in the azimuthal direction do not 
significantly differ from results with higher azimuthal resolution in the nonlinear stage. 

4. Summary and Discussion 

We carried out 3D MHD simulations of the galactic gaseous disk. The numerical results 
indicate that magnetic fields are amplified by MRI. The amplified azimuthal magnetic 
fluxes buoyantly escape from the disk by Parker instability. Azimuthal component of the 
mean magnetic field reverses its direction quasi-periodically. The timescale of the reversal 
is about 300to which corresponds to about 1.5 Gyr. 

The numerical results also show that galactic gaseous disk becomes turbulent and that 
averaged plasma stays at around /3 ~ 10. Plasma /3 becomes lower near the disk surface 
due to the density stratification. Since rotational speed decreases from the disk to the halo 
(see figure Eb), turbulent magnetic fields are stretched around the disk surface, forming 
ordered fields. When the magnetic pressure of the ordered fields becomes as large as (3 ~ 5, 
magnetic flux escapes from the disk due to the buoyancy. 

Figure M schematically illustrates the mechanism of the galactic dynamo obtained from 
numerical simulations. Let us consider a differentially rotating disk threaded by weak seed 
fields, which rotates counter clockwise. As MRI grows in the gaseous disk, positive and 
negative fields whose strengths are comparable are amplified. When the averaged magnetic 
pressure becomes about 10% of the gas pressure in the disk, local magnetic pressure near 
the surface exceeds 20% of the gas pressure and Parker instability creates magnetic loops 
near the disk surface. Since the magnetic field strength parallel to the seed field grows faster 
than that of the opposite field, the counter clockwise fields satisfy the condition of Parker 
instability earlier and buoyantly escape from the disk. Subsequently, clockwise magnetic 
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fields which remain inside the disk are amplified by MRI, and mean azimuthal magnetic 
fields inside the disk reverses their direction. The reversed fields become a seed field for the 
next cycle. 

Similar quasi-periodic reversal of azimut hal fields have b een reported in local 3D MHD 



simulations of a diff erentially rotating disk (Shi et al. 



so long as (5 > 0.1 (IJohansen &: Levin 



20101 ). Although MRI can grow 



20081 ). the saturation level of plasma is about 



~ 1 — 10 when the simulations include the effect of the vertical stratification. The growth 
of MRI saturates when the growth rate of the Parker instability is largest (/3 ~ 1 — 10). It 
indicates that the saturation levels of the magnetic energy inside the disk is determined by 
the Parker instability around the surface layer. 

Our finding about the periodical reversal of the direction of the halo azimuthal field 
would impact on the obse rvation of all sky RM map . The r e are a number of reports ab out 



all sky RM distribu tions ( ITaylor. Stil. &: Sunstrum 



2009 



Braun. Heald fc Beck 



2010 



Han fc Zhang 



20071 ) . Study of all sky RM distribution is essential not only to understand 



the global structure of galactic magnetic fields but also to explore the intergalactic magnetic 



fields in the intergalactic medium 



instance. 



Taylor. Stil. fc Sunstrum 



Akahori & Ryu 2010, 2011, references therein). For 
( 120091 ) suggested that the existence of the halo poloidal 



field is consistent with a dipole field rather than a quadrupole field. 

We calculated the distribution of RM obtained from numerical results (figure ITU]) . The 
position of the observer is assumed to be at r = 8kpc at t — lOOOto- The direction of the 
rotational velocity and azimuthal magnetic fields is counter-clockwise in our numerical 
simulation. Since the angular velocity observed in Milky way is clockwise, we calculated 
the RM distribution by reversing the z-direction (i.e. north is —z direction). From the 
characteristics of MHD, if the initial magnetic field direction were reversed, we can obtain 
the same result except for the direction of magnetic fields. Thus we reversed the color of 
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the RM from the original one. The scale of RM is arbitrarily in this paper. We found 
that magnetic loops produced by Parker instability formed multiple reversals of the the 
sign of RMs along the latitude. These reversals reveal that magnetic fields with opposite 
polarity emerge from the disk periodically. The snapshot shown in figure [10] indicates that 
the distribution of RM is point symmetric with respect to the galactic center, and that this 
feature remains in timescale of ~ 1.5 Gyr (figure 2a). 

Braun et al. (2010) observed RM distributions in nearby galaxies. They also pointed 
out that the field topology in the upper halo of galaxies is a mixture of axi-symmetric spiral 
quadrupole field in thick disks and raidally directed dipole field in halos. They suggested 
that the origin of the dipole components might be a bipolar outflow. Mixed structure 
was also suggested by recent studies of all-sky RM map (see, e.g., figure 8 of Pshirkov 
et al. 2011). The RM distribution in the halo region of our numerical results shows an 
anti-symmetric distribution corresponding to the dipole field. 

Outflows are produced by emerging magnetic fluxs driven by Parker instability. The 
time interval of the buoyant escape of the magnetic flux by Parker instability is about 10 
rotational periods which corresponds to the growth timescale of MRI. Since buoyant escape 
of magnetic flux takes place about 10 cycles during the simulation for 4Gyr at zu/voq = 2, 
the disk becomes turbulent enough, and the influence of the initial field directions become 
negligible. Even when we start simulations with magnetic fields symmetric to the equatorial 
plane, the direction of the halo fields in the upper and lower hemisphere becomes opposite 
at uj/uj = 2 as shown in figure IHK On the other hand, the effect of the initial field topology 
still remains in the outer disk region. Since RM is the sum of the line of sight value, the 
high latitude region (i.e., magnetic loops in nearby disks) still memorizes initial fields. Since 
the outflow velocity may be faster when we include th e effect of the super nova explosions 



( I Wada &: Norman 



20011 ) and the cosmic ray pressure (IKuwabara. Nakamura. fc Ko 



2004 
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Hanasz. Woltanski. fe Kowalik 1120091 ). outer disks may have enough time to create halo 
fields by the buoyant escape of turbulent disk fields, and eliminate the effect of the initial 
fields. 

Our numerical calculation showed that initially non-rotating halo rotates after several 
rotation periods, because the angular momentum was supplied from the disk to halo by 
rising magnetic loops formed by the Parker instability. In fact, HI observation of late type 
spiral galaxy, NGC 6503 indicated that there were two components of HI disk, thin dense 
disk and t hick low density disk , and that the thick low-density disk rotates slower than the 



thin disk ( IGreisen et al 



slower than that of HI ( jBottema 



2009). 



he ro tation speed of the ionized gas in the thick disk is 



19891 ). 



In this paper, we assumed that the interstellar gas is one com ponent gas with 



Nishikori et al. 



tempe rature 2 x 10 5 K in order to compare the result with that by 
( 120061 ). However, in realistic galaxy, the interstellar gas has multi-temperature, multi-phase 
components. In future works, we would like to include the energy input by supernovae, 
multi-temperature structure of the interstellar gas, and cosmic-rays. 
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Table 1: Units of Physical quantities. 



Physical quantity 


Symbol 


Numerical unit 


Length 


w 


1 kpc 


Velocity 


v 


207 km/s 


Time 


to 


4.8 xlO 6 yr 


Density 


Po 


1.6 xl(T 24 g/ cm3 


Temperature 


T 


5.15 xlO 6 K 


Magnetic field 


B 


26 n G 
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Time (t/to) Time (t/t ) 



Fig. 1. — (a) Time evolution of the mass accretion rate at zu/zxj = 2.5. Black curve shows 
model ASYM and gray curve shows model PSYM. (b) Time evolution of the plasma (5 
averaged in the region where 2 < vj/-ujq < 5, \z\/vdq < 1, and < tp < 2n. (c) Time 
evolution of the spatially averaged magnetic energy of model ASYM. Black, dark gray and 
gray denotes azimuthal component, radial component and vertical component, respectively, 
(d) Time evolution of the spatially averaged magnetic energy of model PSYM. Colors are 
same as in (c). The averaged region of the magnetic energy is same as (b). 
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Fig. 2. — Evolutions of the mean azimuthal magnetic fields averaged in the region where 
5 < zu/zu < 6, < (p < 2tt, and < z/zu < 1 (black), and in the region where 5 < 
w/vjq < 6, < ip < 2n, and 1 < z/vjq < 3 (gray), (a) and (b) show time evolution and 
radial distribution, respectively. 
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Fig. 3. — (a) Radial distribution of the specific angular momentum (I = wv^/ {wqVq)) at 
equatorial plane. Black curve shows the distribution at t — lOOOto; and gray curve shows the 
initial condition. Light gray dashed line shows the position of the initial density maximum, 
(b) Spatial distribution of the azimuthally averaged rotation speed at t — lOOOto- 



Fig. 4. — Isosurface of density at t — lOOOio- The region |x| / -Cc7 , \y\/vjQ < 18.5, and \z\/voq < 
7 is plotted. Curves show the mean magnetic field lines. Color denotes direction of the 
azimuthal magnetic fields. Blue-white-red corresponds to negative-positive. 
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Fig. 5. — (a) A typical magnetic loop at t — 1000t - Black circle denotes the starting point 
of integration of the magnetic field line. The distribution of the physical quantities along 
the magnetic field line, (b) vertical velocity, (c) plasma j3, and (d) density, respectively. 
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Time (t/to) Time (t/to) 

Fig. 6. — (a) Time evolution of the azimuthal component of magnetic field averaged in 
the azimuthal direction. White curve shows the location where plasma (3 = 5 above the 
equatorial plane {z > 0). (b) Correlation of the direction of B v below the equatorial plane 
and above it. White shows the region of positive correlation where the direction of azimuthal 
magnetic fields above and below the disk is identical, and black is the region with negative 
correlation. 
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Fig. 7. — (a) w — z distribution of the azimuthal magnetic fields averaged in the azimuthal 
direction at t/to = 1000. (b) w — tp distribution of azimuthal magnetic fields averaged in the 
vertical direction where 1 < z/wq < 1.5. (c) w — tp distribution of the azimuthal magnetic 
fields averaged in < z/zu < 0.5. 
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Fig. 8. — (a) Time evolution of the magnetic energy averaged in the region where 2 < 
w/vjq < 5, \z\/vdq < 1, and < ip < 2n . (b) Time evolution of the azimuthal magnetic field 
averaged in the region where 5 < w/w < 6, < z/vdq < 1, and < (p < 2n. Black, gray, 
and dashed curves show the results for Model ASYM, ASYM64, and ASYM256, respectively. 



- 29 - 

(a) 




Fig. 9. — Schematic drawing of the mechanism of MHD Dynamo. 
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Fig. 10- 



Rotation measure distribution obtained from numerical results. 



